REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-01-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  averege  1  hour  per  responso,  including  the  time  for  reviewing  instructions,  searching  existing  dete  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  eny  other  espect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden  to  Deportment  of  Defense,  Washington  Headquarters  Services  Directorate  for  Information  Operetions  and  Reports 
(0704-0188)  1215  Jefferson  Devis  Highwey,  Suite  1204.  Arlington  VA  22202-4302.  Respondents  should  be  ewere  that  notwithstanding  any  other  provision  of  lew.  no  person  shell  be 
subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  e  currently  valid  OMB  control  number 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

15-11-2008 

REPRINT 

4.  TITLE  AND  SUBTITLE 


On  the  anticritical  temperature  for  spacecraft  charging 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 

62601F 


6.  AUTHORS 

Shu  T.  Lai  and  Maurice  Tautz* 


5d.  PROJECT  NUMBER 

5021 


5e.  TASK  NUMBER 

RS 


5f.  WORK  UNIT  NUMBER 

A1 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  RVBXT 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-RV-HA-TR-2008-1 142 


9  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/RVBXT 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  distribution  unlimited. 


1 3.  SUPPLEMENTARY  NOTES 

Reprinted  from  Journal  of  Geophysical  Research.  Vol  113,  A 1 121 1,  doi :  1 0  1 029/2008J A0 13161,  2008  ©  2008  American  Geophysical  Union 
*AER/Radex,  Inc.,  Lexington,  MA 


14.  ABSTRACT  In  recent  years,  evidence  has  been  found  for  the  existence  of  a  critical  temperature  for  the  onset  if  spacecraft  charging  to 

high  voltages  High-voltage  charging  affects  scientific  instruments  on  board  and  is  related  to  spacecraft  anomalies.  However,  less  attention  has 
been  given  to  low-voltage  charging,  which  can  also  affect  scientific  experiments  on  board  and  is  relevant  to  surface  chemistry.  There  also  can 
exist  an  anticritical  temperature  for  low-voltage  spacecraft  surface  charging.  Ambient  electrons  at  very  low  tremperatures  tend  to  cause  negative 
surface  charging,  albeit  at  low  voltages,  and  as  the  electron  temperature  increases,  the  charging  ceases  at  a  critical  value  depending  on  the  surface 
material  We  present  the  theory  and  numerical  results  of  anticritical  temperatures  for  typical  surface  materials  in  Maxwellian  space  plasmas.  The 
change  in  anticritical  temperature  due  to  a  low-incident-energy  enhancement  of  the  electron  backscatter  yield,  consistent  with  recent  measure¬ 
ments,  is  discussed.  Approximate  expressions  for  the  anticritical  temperature  upper  limits  are  given  on  the  basis  of  Taylor  expansions  at  low 
temperature  of  the  charging  onset  equation.  It  is  shown  that  the  existence  of  the  anticritical  temperature  slightly  modifies  the  possible  triple-root 
configurations  in  the  flux-voltage  characteristic  curve  for  a  material  The  surface  charging  effect  of  a  Maxwellian  plasma  with  flux  components 
spanning  the  anticritical  and  critical  temperatures  is  considered.  A  comparison  with  an  empirical  low-voltage  charging  curve  is  given 


15.  SUBJECT  TERMS 

Spacecraft  charging 
Secondary  electrons 


Critical  temperature 
Geosynchronous  satellite 


Anticritical  temperature 


Space  plasmas 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

15 

Shu  T.  Lai  In  recent  years,  evidence  has  heen 

UNCL 

UNCL 

UNCL 

Uni 

19B.  TELEPHONE  NUMBER  (Include  area  code) 

Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std  239  18 


AFRL-RV-HA-TR-2008-1 142 

JOURNAL  OF  GEOPHYSICAL  RESEARCH,  VOL.  113,  A1 1211.  doi:  1 0. 1 029/2008 JAO 13161,  2008 


anticritical  temperature  for  spacecraft  charging 

Shu  T.  Lai'  and  Maurice  Tautz2-3 

Received  1 1  March  2008:  revised  5  June  2008;  accepted  18  August  2008;  published  15  November  2008 

[t]  In  recent  years,  evidence  has  been  found  for  the  existence  of  a  critical  temperature 
for  the  onset  of  spacecraft  charging  to  high  voltages.  High-voltage  charging  affects 
scientific  instruments  on  board  and  is  related  to  spacecraft  anomalies.  However,  less 
attention  has  been  given  to  low-voltage  charging  which  can  also  affect  scientific 
experiments  on  board  and  is  relevant  to  surface  chemistry.  There  also  can  exist  an 
anticritical  temperature  for  low-voltage  spacecraft  surface  charging.  Ambient  electrons  at 
very  low  temperatures  tend  to  cause  negative  surface  charging,  albeit  at  low  voltages, 
and  as  the  electron  temperature  increases,  the  charging  ceases  at  a  critical  value  depending 
on  the  surface  material.  We  present  the  theory  and  numerical  results  of  anticritical 
temperatures  for  typical  surface  materials  in  Maxwellian  space  plasmas.  The  change  in 
anticritical  temperature  due  to  a  low-incident-energy  enhancement  of  the  electron 
backscatter  yield,  consistent  with  recent  measurements,  is  discussed.  Approximate 
expressions  for  the  anticritical  temperature  upper  limits  are  given  on  the  basis  of  Taylor 
expansions  at  low  temperature  of  the  charging  onset  equation.  It  is  shown  that  that  the 
existence  of  the  anticritical  temperature  slightly  modifies  the  possible  triple-root 
configurations  in  the  flux-voltage  characteristic  curve  for  a  material.  The  surface  charging 
effect  of  a  Maxwellian  plasma  with  flux  components  spanning  the  anticritical  and 
critical  temperatures  is  considered.  A  comparison  with  an  empirical  low-voltage  charging 
curve  ts  given. 

Citation:  Lai,  S.  T.,  and  M.  Tautz  (2(X)8),  On  the  anlicriiical  lempcralure  lor  spacecraft  charging,  J.  Geophys.  Res..  1 1 3.  A 1 1 2 1 1 , 
doi:  10. 1029/2008 JAO  13 161 . 


Click 

Here 

Full 

Article 
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1.  Introduction 

[2]  In  recent  years,  evidence  has  been  found  in  support  of 
a  entieal  temperature  F*  of  ambient  electrons  for  the  onset 
of  spacecraft  charging  at  geosynchronous  altitudes.  The 
critical  temperature  is  the  threshold  level  for  negative 
voltage  charging,  based  solely  on  electron  currents,  assum¬ 
ing  a  Maxwellian  distribution  of  incident  electrons,  with 
accompanying  secondary  and  baekseatter  interactions  in  the 
satellite  surface  material.  The  net  current,  integrated  over 
energy,  depends  on  the  temperature  of  the  Maxwellian 
distribution.  The  threshold  temperature  oeeurs  when  there 
is  a  balance  between  the  primary  electron  currents  entering 
the  material  and  the  seeondary'baekseatter  electrons  escap¬ 
ing  from  it.  The  solution  to  this  condition  yields  F*  the 
critical  temperature  for  the  onset  of  spacecraft  charging  [Lai 
et  al. ,  1982,  1983;  Laframboise  et  ai ,  1982;  Laframboise 
and  Kamitsuma ,  1983;  Lai ,  1991a,  1991b;  Lai  and  Della- 
Rose ,  2001;  Lai  and  Tautz ,  2006] 

[3]  For  many  materials,  because  of  the  shape  of  the 
secondary  and  baekseatter  emission  curves,  there  is  a  second 
solution  at  a  lower- temperature  TA.  The  solution,  TA ,  has  the 
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property  that  surfaee  charging  to  negative  voltages  oeeurs 
below  Ta  but  is  suppressed  above  it  Because  this  effect  is 
the  opposite  to  that  of  the  critical  temperature,  we  call  TA  the 
anticritical  temperature.  In  sections  2  and  3,  we  outline  the 
theory  of  the  entieal  and  anticritical  temperatures.  Approx¬ 
imate  calculations  of  the  anticritical  temperature  upper  limits 
are  presented  in  seetion  4.  In  section  5  we  outline  the  physics 
of  triple  root  formation  and  in  sections  6  and  7  we  study  the 
problem  of  triple-root  jumps  and  equilibrium  spacecraft 
potentials  in  a  double  and  triple  Maxwellian  plasma,  assum¬ 
ing  that  both  Ta  and  F*  exist  for  a  given  surfaee  material. 
Seetion  8  gives  a  discussion  of  the  charging  effect  of  a  tnple- 
Maxwellian  spaee  plasma,  having  temperature  components 
spanning  the  anticritical  and  critical  temperatures,  with 
emphasis  on  the  low-voltage  spacecraft  charging  region.  In 
seetion  9  we  give  a  conclusion. 

2,  Existence  of  the  Critical  Temperature 

[4]  We  first  outline  the  concept  of  critical  temperature.  The 
ambient  electron  eurrent  often  exceeds  that  of  ambient  ions  by 
nearly  2  orders  of  magnitude  [Reagan  et  ai ,  1983]  because  of  the 
mass  difference  and  thus  the  ion  contribution  can  be  neglected. 
The  onset  of  spacecraft  charging  at  equilibrium  is  determined  by 
the  balance  of  incoming  and  outgoing  electron  fluxes 

J  </’vv/(v)  =  J  Jivvn[6(E)  +  »/(£)]F(v),  (la) 
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where  E(v)  is  the  ambient  electron  velocity  distribution,  v„ 
is  the  normal  component  of  the  electron  velocity,  and  both 
sides  are  actually  definite  integrals.  One  can  write  vn  = 
v’cosfl  in  spherical  coordinates  (v,  6 ,  0),  where  0  is  from  0  to 
7t/2  and  &  from  0  to  27t.  Writing  the  integrals  explicitly, 
equation  (la)  is  of  the  form 

r-K  rlffl  r2n 

/  dvv1  cIO  sin9  /  d(pvnF(v) 

Jo  Jo  Jo 

/•x  r^/2  r  In 

=  /  dvv1  /  dOsinO  /  d(pvn[6(E)  +  r/(£)]E(v).  (lb) 
Jo  Jo  Jo 

In  equation  (la),  the  secondary  and  baekseattered  electron 
yields,  6(E)  and  r/(E)  respectively,  are  measured  in  the 
laboratory  as  functions  of  electron  energy,  E,  where  E  ~ 
(1/2  )mv2. 

[5]  To  solve  equation  (la)  for  a  given  surface  material, 
one  needs  to  input  the  functions  6(E)  and  77(E).  It  is 
convenient  to  use  the  fitted  /i(E)  function  obtained  by 
Sanders  and  fnouye  [1978] 

6(E)  =  c[exp(-E/a)  -  exp(  E/b)}  (2) 

where  a  =  4.3  £max^  b  ~  0.367 Emax,  and  c  =  1.37  6max. 

[6]  For  the  baekseattered  electron  yield  r](E ),  Prokopenko 
and  Laframboise  [1980]  gave  the  fitted  function 

r](E)  =A  -Bexp(-CE),  (3) 

where  the  parameters  Emax,  6 max,  Ay  B ,  and  C  depend  on  the 
surfaee  material. 

[7]  Converting  to  spherical  coordinates  (E,  0 ,  p)y  one  can 
integrate  equation  (la)  over  the  energy  and  the  two  angular 
variables.  For  normal  incidence,  the  angles  on  both  sides  of 
the  equation  eaneel  out.  After  some  algebra,  one  obtains  the 
following  equation: 


[  dEEf(E)  [  dEE[6(E)  +  r/(E)l/*(E),  (4) 

Jo  Jo 

where  /(E),  obtained  from  E(v)  is  the  electron  distribution 
with  the  velocity  expressed  in  terms  of  energy.  For 
Maxwellian  distribution,  j(E)  is  given  by 

f{E)  =  (rt/2mkT)i/2  cxp(- E/kT),  (5) 

where  n  is  the  electron  density,  m  the  electron  mass,  T  the 
electron  temperature,  and  A  the  Boltzmann  constant. 
(Strictly  speaking,  equation  (5)  is  not  an  energy  distribution. 
There  is  an  energy  distribution,  which  is  differently  defined 
[see,  e.g.,  A layer  and  Mayer ,  1963],  but  we  will  not  need  to 
use  it  here.)  Sinee  the  plasma  density  n  is  multiplicative,  it 
cancels  on  both  sides  of  equation  (4)  Therefore,  the 
threshold  condition  is  independent  of  the  electron  density  n 
and  depends  on  the  temperature  T  only.  The  condition  (4) 
ean  be  written  in  the  equivalent  form 

JT  dEE{6(E)  +  ii(E)\f(E)  , 

— Tmd) - '•  <6,) 


In  shorthand  notations,  equation  (6a)  can  be  written  as 

{6  +  v)  =  1,  (6b) 

where  the  notations,  (  and  ),  denote  averaging  in  the  sense 
of  equation  (6a). 

[8]  Substituting  equations  (2)  and  (3)  into  equation  (4), 
with  the  distribution  function  /(E)  given  by  equation  (5)  and 
integrating  over  E,  one  obtains  a  sixth-order  algebraic 
equation  for  AT,  corresponding  to  equation  (4) 

c[(l  +  kT/a)  2-(l  +kT/b)~2}  +  A  -  B(C  kT  +  1)'2=  1  (7) 

A  solution  of  equation  (7)  can  be  used  to  specify  the  critical 
temperature  T*  [ Lai  et  al. ,  1 982]  for  the  onset  of  spacecraft  charging. 

[9]  For  electrons  coming  in  at  various  angles,  one  needs 
to  use  angle-dependent  6  and  rj  functions.  Angle-dependent 
6  and  rj  functions  have  been  given  by  Darlington  and 
Coslett  [1972]  and  Prokopenko  and  Laframboise  [1980], 
respectively.  The  algebra  becomes  more  complicated  but  the 
physics  remains  unchanged.  Further  details  on  the  critical 
temperature  arc  given  by,  for  example,  Lai  and  Della-Rose 
[2001]  and  Lai  and  Tautz  [2006]  and  will  not  be  repeated 
here.  For  comparison  purposes,  we  include  Table  1  listing 
the  critical  temperature  T*  calculated  for  several  typical 
surface  materials,  based  on  material  properties  given  by 
Laframboise  and  Kamitsuma  [1983].  One  can  see  that  kT* 
is  typieally  in  the  few  keV  range. 

[10]  There  are  recent  advances  in  secondary  electron  yield 
formulae  derived  by  using  laboratory  measurements  [e.g., 
Katz  et  al, ,  1986;  Cazaux ,  2001,  2006;  Lin  and  Joy ,  2005; 
J.  R.  Dennison,  Physics  Department,  Utah  State  University, 
manusenpt  in  preparation].  We  are  doing  a  comparative 
study  to  assess  the  effects  of  the  various  formulae  on  the 
critical  temperature.  Whenever  an  improved  formula  of  6(E) 
or  77(E)  becomes  available,  one  can  substitute  it  in  equation 
(4)  for  updating  the  numerical  value  of  the  entieal  temper¬ 
ature  and  the  anticntieal  temperature. 

3.  Existence  of  the  Anticritical  Temperature 

[11]  Figure  1  shows  the  yields  for  the  sample  material 
gold  as  a  function  of  incident  energy  in  the  energy  range  0  < 
E  <  5  (keV).  The  lower  curve  gives  the  baekseatter  yield, 


Tabic  1.  Critical  Temperature  7*a 


Material 

Isotropic 

Normal 

Mg 

0.4 

Al 

0.6 

- 

Kapton 

0.8 

0.5 

Al  Oxide 

2.0 

1.2 

Teflon 

2.1 

1.4 

Cu-Bc 

2.1 

1.4 

Glass 

2.2 

1.4 

Si02 

2.6 

1.7 

Silver 

2.7 

1.2 

Mg  Oxide 

3.6 

2.5 

Indium  Oxide 

3.6 

2.5 

Gold 

4.9 

2.9 

Cu-Bc  (Activated) 

5.3 

3.7 

MgF2 

10.9 

7.8 

aUnits  in  keV  for  kT*. 
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ELECTRON  ENERGY  E(KeV) 

Figure  1.  Yield  versus  energy  for  gold. 


the  middle  curve  the  secondary  emission  and  the  top  curve 
the  summed  yield,  b(E)  +  r](E).  It  can  be  seen  that  at  high 
energies,  the  total  yield  goes  below  unity,  giving  a  negative 
current  contribution,  and  the  integration  over  energies  can 
produce  current  balance  (see  equation  (4))  at  the  critical 
value  7*.  It  is  apparent  from  Figure  1  that  the  “average” 
yield  (6  +  77)  integrated  over  energy  (equation  (3))  can  also 
go  below  unity  at  low  energies  and  that  there  could  be  a 
second  solution.  We  shall  call  the  unity  crossing  point  ({b  + 
rj)  -  1 )  at  low  temperature  TA ,  the  antieritical  temperature. 
The  existence  and  value  of  TA  depend  critically  on  the  low- 
energy  behavior  of  the  yield  curves. 

[12]  Table  2  lists  the  numerically  calculated  values  of 
antieritical  temperature  TA  for  some  typical  surface  materi¬ 
als,  based  on  the  parameter  coefficients  of  Sanders  and 
Inouye  [1978]  and  Prokopenko  and  Laframboise  [1980].  It 
is  observed  that  the  AT  1  values  arc  mostly  below  20  eV  for 


normal  incidence  and  below  10  cV  or  nearly  zero  for 
isotropic  incidence  Such  low  values  of  antieritical  temper¬ 
ature  are  easily  overlooked.  From  Table  2,  we  sec  that  the 
b(E)  coefficient  of  some  materials,  such  as  gold  and  silver, 
has  two  properties:  (1)  high  values  of  Emax  and  (2)  near 
unity  value  of  ^max.  These  two  properties  give  higher  values 
of  TA.  The  maximum  value  found  is  288  cV  for  silver  at 
normal  incidence, 

[13]  Since  the  secondary  and  backscattcred  electron  emis¬ 
sion  coefficients  b(E)  and  rfE)  arc  obtained  by  fitting  data 
obtained  in  the  laboratory,  the  signal-to-noise  ratio  often 
decreases  as  the  energy  E  approaches  zero,  rendering 
inaccuracies  to  the  fitted  data  in  the  low-energy  range 
Therefore,  the  coefficients  b(E)  and  tj(E)  should  not  be 
regarded  as  highly  accurate  for  small  E  In  Table  2,  anti- 
critical  temperatures  kTA  below  1  eV  are  flagged  and  no 
value  is  given.  It  is  also  possible  that  the  secondary  and 
baekseatter  emissions  are  so  small  that  (b{E)  +  J]{E)}  is  less 
than  one  at  all  energies.  In  that  case,  there  would  not  be  a 
solution  for  either  the  critical  or  antieritical  temperature  (see 
materials  Mg  and  A1  in  the  Tables  1  and  2,  respectively)  and 
the  material  should  always  charge. 

[14]  The  existence  of  the  antieritical  temperature  TA  stems 
from  the  properties  of  the  secondary  electron  emission 
coefficient  b(E)  and  baekseatter  coefficient  ifE).  The  inci¬ 
dent  electrons  must  have  enough  energy  to  generate  a 
secondary  electron  At  sufficiently  low-energy  E  of  the 
incoming  electrons,  not  many  exited  electrons  are  ereated 
in  the  material.  The  result  is  that  the  number  of  incoming 
(primary)  electrons  exceeds  that  of  the  outgoing  secondary 
electrons  and  the  secondary  emission  yield  goes  to  zero  as  E 
approaches  zero.  This  is  not  necessarily  true  for  the  baek- 
scattcrcd  electrons.  In  Figure  1,  the  baekseatter  yield  goes  to 
a  small  positive  value,  which  is  set  by  a  fit  to  extrapolated 
measured  data  [Prokopenko  and  Laframboise ,  1980]  In 
terms  of  the  baekseatter  parameters,  the  yield  at  E  =  0  is  A 

B ,  which  for  gold  is  0.124.  However,  there  arc  experimental 
indications  that  low-energy  backscattcred  particles  have  a 
reflection  coefficient  of  one  [Cirnino  et  ai ,  2004]  as  E  goes 
to  zero.  This  result  can  also  be  derived  using  a  quantum 
mechanics  model  considering  a  plane  wave  electron  inci- 


Table  2.  Antieritical  Temperature  T4  for  Typical  Surface  Materials3 


Material 

f 

^max 

^max 

A 

B 

C 

lsotropicb 

Normal6 

Mg 

0.25 

0.92 

0.1460 

0.0250 

0.3440 

_ 

_ 

Al 

0.30 

0.97 

0.1568 

0.0303 

0.3431 

— 

— 

Kapton 

0.15 

2.10 

0.0700 

0  0000 

0.0000 

£ 

0.014 

AUO, 

0.30 

2.60 

0.1238 

0.0172 

0.3455 

£ 

0.020 

Teflon 

0.30 

3.00 

0.0900 

0  0000 

0.0000 

e 

0.017 

Cu-Bc 

0.30 

2.20 

0.3136 

0.0692 

0.6207 

£■ 

0.020 

Glass 

0.35 

2.35 

0.2000 

0.0420 

0.4100 

£ 

0.015 

SiO, 

0.42 

2.50 

0.1238 

0.0172 

0.3435 

0.024 

0.029 

Silver 

0.80 

1.00 

0.3900 

0.2890 

0.6320 

0.184 

0.288 

Magnesium  Oxide 

0.40 

4.00 

0.1238 

00172 

0.3435 

e 

0.016 

Indium  Oxide 

0.80 

1.80 

0.2750 

0.0250 

0.5400 

0.077 

0.108 

SCATHA  Boomai 

0.59 

1.86 

0.4380 

0.3250 

0.6130 

0.041 

0.059 

Gold 

0.80 

1.45 

0.4802 

0.3566 

0.6103 

0.096 

0.114 

Cu-Bc  (Activated) 

0.40 

5.00 

0.3)36 

0.0692 

0.6207 

C 

0.010 

MgF, 

0.85 

6.38 

0.1238 

0.0172 

0.3435 

c* 

0.019 

“Antieritical  temperature  TA  given  in  kcV. 

6Hcrc  means  that  there  is  no  TA ,  and  e  indicates  that  TA  is  below  1  cV. 
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ELECTRON  TEMPERATURE  *T(keV) 


Figure  2a.  Current  balance  equation  (b  +  rf)  =  1  for  the 
antieritieal  temperature  TA  with  and  without  enhanced 
baekseatter  for  the  surface  material  gold.  We  have  used 
£*0  -  0.05  keV  for  the  enhancement  fall-off  parameter. 

dent  on  a  negative  potential  step  [see,  e.g.,  Bransden  and 
Joachain ,  1089]. 

[15]  We  have  calculated  the  yields  for  gold  in  the  low- 
energy  range  0  <  E  <  0.3  (keV),  with  an  added  low-energy 
baekseatter  enhancement.  The  material  properties  are  the 
same  as  in  Figure  1 ,  except  for  the  additional  term 

A»/=(l  -*+fi)«xp(- (8) 

where  parameter  £()  specifies  the  enhancement  fall-off  rate. 
With  this  change,  r](E)  +  Ar)(E)  — ►  1  as  E  — » ►  0  and  the  same 
asymptotie  baekseatter  yield  is  obtained  at  large  E.  Experi¬ 
mental  data  indicates  that  £0  is  about  0.05  1 .5  keV  for  gold 
[Jablonski  et  ai ,  1989;  Jablonski  and  Jiricek ,  1996]  and 
0.05  keV  has  been  used.  The  effect  of  the  added  term  is  to 
produce  a  minimum  in  the  summed  yield  eurve.  If  this 
minimum  is  less  than  1 .0,  an  antieritieal  temperature  TA  exists. 

[16]  We  show  in  Figure  2a  the  energy  averaged  coeffi¬ 
cient  a  =  (6  +  rj  +  Arf)  for  gold  as  a  function  of  kT  (keV), 
with  (solid  line)  and  without  (dashed  line)  the  additional 
baekseatter  enhancement  term,  (A rf) 

(A,,)  =  (\-A  +  B)/(\+kT/E0)2.  (9) 

For  reference,  the  primary  eurve  corresponding  to  the  no 
baekseatter  ease  is  shown  as  a  dash  line.  The  antieritieal 
temperature  is  determined  by  the  intersection  of  the  curves 
with  the  horizontal  “negative-to-positive”  transition  line  at 
<7  =  1.0.  With  the  added  term  TA  was  determined  numerically 
to  be  0.0915,  which  can  be  compared  with  0.1143  for  the 
unmodified  ease.  The  TA  value  is  shifted  to  lower  kT, 
because  the  <7  eurve  is  elevated  by  the  baekseatter  so  that  the 
intersection  with  the  horizontal  test  line  oeeurs  earlier.  The 
values  of  antieritieal  temperature  given  in  Table  2  thus 
represent  upper  limits. 


[17]  Near  the  critical  temperature,  there  would  also  be  a 
shift,  in  this  ease  to  higher  kT,  though  it  would  be  small 
since  A rj(E)  falls  off  rapidly  with  £.  For  the  ease  of  gold, 
we  obtain  2.913  for  the  shifted  critical  temperature  and 
2.911  for  the  unmodified  ease. 

[is]  For  some  materials,  the  baekseatter  enhancement 
could  be  strong  enough  so  that  there  is  no  zero  crossing 
in  the  net  averaged  flux,  (1  6  rf),  eurve  and  no 

antieritieal  temperature  would  occur.  We  do  not  ealeulate 
a  table  of  modified  antieritieal  temperatures  because  the 
parameter  £0  is  not  generally  known  experimentally.  As 
indicated  in  Figure  2b,  when  £0  goes  to  zero,  we  recover  the 
unmodified  antieritieal  temperature.  As  £(>  increases,  the 
antieritieal  temperature  decreases  until  at  some  point  (here 
£0  ~  0. 1 3)  it  disappears.  This  type  of  behavior  would  oeeur 
with  other  materials  and  the  magnitude  of  the  temperature 
shift  £  would  depend  on  the  specific  material  parameters. 
The  shift  can  be  estimated  from  the  current  balance  equa¬ 
tion,  if  we  expand  about  kT  of  the  no  enhancement  ease.  For 
example,  we  get  for  gold 

A Ta  =  -Ar)/(da/dkT)  «  -0.081/3.4  =  0  024.  (10) 

whieh  agrees  approximately  with  the  antieritieal  tempera¬ 
ture  difference  (-0.023)  found  numerically.  At  high  kT,  in 
the  critical  temperature  region,  the  derivative  doklkT 
changes  sign  and  the  enhancement  term  is  small,  leading 
to  a  slight  positive  shift. 

[19]  There  are  other  baekscattenng  formulae  and  meas¬ 
urements  published  in  the  literature.  They  ean  slightly  affect 
the  numerical  results.  Recently,  Cimino  et  ai  [2004]  mea¬ 
sured  the  baekscattcring  of  gold  and  found  deviation  at  near 
zero  energy.  In  view  of  this  uncertainty,  we  have  done  a 
study  using  the  baekscattenng  coefficient  for  gold  as  a 
parameter  going  over  the  expected  range  from  0. 1 24  to  1 
at  near  zero  energy.  Since  no  measurement  of  other  maten- 
als  is  available  at  near  zero  energy  at  this  time,  we  did  not 
pursue  this  aspect.  The  baekscattcring  equation  does  not 
change  the  concept  of  our  paper.  At  this  time,  Cimino’s 


Figure  2b.  Antieritieal  temperature  TA  versus  the  baek¬ 
seatter  parameter  £0,  with  gold  as  the  surface  material. 


4  of  15 


A11211 


LAI  AND  TAUTZ:  ANT1CR1TICAL  TEMPERATURE 


A11211 


Table  3.  Anticritical  Temperature  and  Three  Successive  Approximations* 


Material 

krA 

kTcs  ±  % 

kT\  ±  % 

kTy 

±% 

kTi  ±  % 

Cu-Bc  (activated) 

0.0099 

0.0117  +18.7 

0.0101  +2.35 

0.0099  +0.04 

0  0099  +0.09 

Kapton 

0.0141 

0.0105 

25.6 

0.0150  +01 

0.0140 

0.20 

0.0141  +0.64 

Mg  Oxide 

0.0154 

0.0146 

4.8 

0.0159  +3.50 

0.0154  +0.04 

0.0154  +0.16 

Teflon 

0.0168 

0.0146 

12.7 

0.0176  +4.82 

0.0168  +0.01 

0.0168  +0.27 

MgF2 

0.0192 

0.0195  +1.6 

0.0196  +2.16 

0.0192  +0.05 

0.0192  +0.08 

ai2o, 

00197 

0  0169 

14.3 

0.0208  +5.48 

00197 

0.02 

0.0198  +0.35 

Cu-Bc 

0.0197 

0  0200  +1.6 

0.0207  +5.53 

0.0196 

0,02 

0.0197  +0.34 

Glass 

0.0243 

0.0218 

10.2 

0.0257  +5.74 

0.0243 

0.03 

0.0244  +0.39 

Si02 

0.0291 

0.0246 

15.4 

0.0307  +5.72 

0.0291 

0.03 

0.0292  +0.39 

SCATHA  Boomat 

0.0595 

0.0465 

21  9 

0.0646  +8.58 

0.0593 

0.27 

0.0598  +0,49 

Indium  Oxide 

0.1082 

0.0837 

22,6 

0.1177  +8.85 

0.1074 

0.73 

0.1094  +1.12 

Gold 

0.1143 

0.0808 

29.3 

0.1281  +12.10 

0.1132 

0.96 

0.1146  +0.31 

Silver 

0.2873 

0.1172 

59.2 

0.3237  +12.66 

0.2615 

8.99 

0.2910  +1.29 

■Units  of  the  temperatures  (multiplied  by  k)  arc  expressed  in  kcV.  Here  %  is  the  percent  difference  from  TA ,  TA  is  the  numerical  solution  with  all  terms,  T0  ±  %  is 
the  analytical  solution  neglecting  the  ion  term,  T\±%  is  the  solution  to  first  order,  T2  ±  %  is  the  solution  to  second  order,  and  7\  ±  %  is  the  solution  to  third  order 


measurement  is  the  only  one  available  at  near  zero  energy. 
One  can  update  the  numerical  values  when  measurements 
for  more  materials  are  available. 

4.  Approximate  Formula  for  Anticritical 
Temperature 

[20]  A  simple  approximate  formula  of  TA  in  Tabic  2  is 
handy  for  practical  purposes.  To  seek  an  approximation,  we 
neglect  the  electron  backscattcring  contribution  compared 
with  that  of  secondary  emission.  With  this  change,  the 
current  balance  equation,  equation  (7),  becomes 

<•[(1  2  -  (I  +kT/b)  2]  =  I  (II) 

By  making  a  Taylor  expansion  in  AT,  and  keeping  only  the 
lowest  terms,  we  get  the  approximate  solution  AF0 

|l2‘> 


or  in  terms  of  Tmax  and  <*>max 


retained,  but  all  powers  of  the  kTib  term  arc  kept  This 
procedure  leads  to  the  cubic  equation  shown  below 

„  _  1  +  2tr(2  _  ,)  +  («  -  i?)  -  m> |  -  0.  ( 13) 


where  we  have  defined  a  -  1  +  (A  B  l  )ic  and  j3  =  1  la 
BC/c.  This  equation  can  be  solved  to  successive  orders  by 
keeping  the  two  lowcst-order  terms,  the  first  three  terms  or  all 
four  terms.  The  solutions  are  denoted  by  kT\,  AT2,  and  AT* 
and  the  results,  along  with  the  percent  errors,  arc  shown  in 
Table  3.  It  is  evident  that  the  linear  estimate  kT\  is  good  to 
about  the  10%  level  and  tends  to  give  high  values.  The 
quadratic  result  AT>  is  accurate  to  about  1%,  but  it  too  has 
10%  errors  in  the  higher-AT  range.  In  order  to  get  the  errors 
down  to  approximately  the  1  %  level  throughout,  one  needs  to 
solve  the  cubic  equation  for  A  7V  However,  if  one  requires  just 
a  practical  “ballpark”  estimate  of  the  anticritical  temperature 
the  linear  result  suffices  and  it  is  explicitly 


b{  1  —  o )  A7o(l  -  A  +  B) 

2(o  b3)  ~  1  +  2kT0[BC  (I  A  f  B)/b] 


(14) 


AT0  =  0.147^.  (12b) 

^max 

This  simple  formula  is  valid  at  small  kT  and  with  no 
backscattcr  component. 

[21]  The  AT0  results  obtained  from  equation  (12b)  arc 
given  in  Table  3  for  a  list  of  available  surface  materials  (sec 
Table  2  for  the  properties).  In  Table  3,  the  first  column 
specifics  the  material.  The  next  column,  labeled  kTA,  is  the 
numerical  result.  The  next  two  columns  show  AT0  and  its 
percent  error.  It  can  be  seen  that,  at  low  anticritical  temper¬ 
atures  (AT  <01  keV)  the  estimate  is  accurate  to  about  20%. 
But  at  higher  temperatures  (AT  >  0.1  kcV),  the  estimate  is 
considerably  worse  (silver  is  ~-60%).  The  main  reason  for 
this  poor  AT0  behavior  can  be  traced  to  the  kTib  term.  Since  b 
itself  is  typically  quite  small,  the  approximation  of  neglecting 
the  higher-order  terms  progressively  fails  as  kT  increases. 

[22]  We  can  obtain  a  better  estimate  of  AT  if  we  develop  a 
hybrid  expansion  of  equation  (7).  By  this  we  mean  that  the 
terms  kT/a  and  CAT  arc  expanded  and  the  first-order  terms 


This  formula  contains  AT0  and  the  extra  terms  lead  to  better 
accuracy  than  equation  (12b).  It  differs  from  the  pure  AT0 
estimate  in  that  backscattcr  coefficients  A,  B,C  are  used  and 
kT/b  has  been  treated  more  fully.  If  E{)  is  known  one  can 
further  estimate  the  effect  of  backscattcr  enhancement  by 
employing  the  shift  equation  (10). 

5.  Triple-Root  Situation  With  Two  Maxwellian 
Components 

[23]  Consider  a  plasma  electron  population  which  is  a 
superposition  of  two  Maxwellian  distributions  w  ith  temper¬ 
atures  T )  and  T2  A  necessary  condition  [Lai,  1991a]  for  a 
triple-root  situation  to  occur  is 

T]  <  T*  <  r2,  (15) 

where  T*  is  the  critical  temperature  Before  we  introduce  a 
new  degree  of  complication,  the  anticritical  temperature  TA , 
into  the  inequality  equation  (15),  let  us  first  review  the 
underlying  idea. 
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[24]  Consider  a  double  Maxwellian  distribution  j(E) 
where  wc  have  expressed  the  velocity  in  terms  of  the  energy 
E.  The  electron  distributions  are 

/(£)=/,(£)  +f2(E).  (16) 


/,(£)  =  n\  (m/2nkTe,\ )3  2  e\p(—E/kTe,\ ),  (17) 


ME)  =  n2{m/2nkTe,2y  2  exp(-E/kTe.2),  (18) 

where  Te\  and  Te2  denote  the  first  and  the  second 
Maxwellian  electron  temperatures  respectively.  Three 
analogous  equations  can  be  written  for  double  Maxwellian 
ions,  with  the  subscript  e  replaced  by  i.  The  temperature  in  a 
Maxwellian  distribution  is  defined  as  the  inverse  of  the 
slope  of  the  graph  of  log  J[E)  as  a  function  of  E.  We  assume 
that  Te  |  <  Te2  and  TiA  <  Tl2.  for  simplicity,  we  assume 
charge  neutrality,  n  =  ne  =  //,.  (In  the  discussion  below, 
variable  densities  arc  indirectly  accounted  for  by  the  varying 
fluxes.)  We  temporarily  neglect  the  ions  for  low  potential 
magnitude.  If  both  temperatures  Te  l  and  Te  l  are  below  T*, 
the  net  flux,  integrated  over  energy  and  including  secondary 
and  backscatter  contributions,  would  be  positive.  The 
surface  would  charge  to  a  positive  potential  and  there  is 
no  competition.  Likewise,  if  both  Te  ]  and  Te  2  are  above  T*, 
the  surface  would  charge  to  a  negative  potential  and  there  is 
again  no  competition  Therefore,  it  is  necessary  to  satisfy 
the  inequality  (15)  for  the  two  distributions  to  compete  with  each 
other.  Such  competition  can  produce  a  triple  root  situation. 

[25]  Let  us  now  write  down  the  full  expressions  of  the 
fluxes  collected  by  a  spacecraft  at  a  negative  voltage  0.  The 
total  flux  JT  is  given  by 

Mq)=M0)+J2(O),  (19) 

where  each  term  is  the  sum  of  the  separate  contributions 
from  the  electrons  and  ions 

J\(4>)=J'.\(4>)+Ji.\{4>)  (20a) 


-4(0)  ='4.2(0)  "b  *4.2  (0)* 

We  have  explicitly  the  terms 

-M*)  =  A.i(0)[l  -  «  + 1)1 '<p(-^) 
+yufl»[l  ~Tf~ 

"'ll 


(20b) 


(21) 


where 


/JC,  x  ^dEEfn(E)[8(E)  +  q(E)\ 

+  = - P  wzr(T?\ - 

Jo  dEEUE) 


(«=  E  2),  (23) 


qe=  -e,  and  </,  =  +e. 

[26]  In  equations  (21)  and  (22 ),yV<w(0)  and  y,„(0)  are  the 
ambient  electron  and  ion  fluxes,  6(E)  and  r](E)  arc  the 
emission  coefficients  for  secondary  electrons  [e.g., 
Sternglass ,  1954a;  Sanders  and  Inouye ,  1978]  and  back- 
scattered  electrons  [e.g.,  Sternglass ,  1954b;  Prokopenko 
and  Laframboise ,  1980],  and  e  the  elementary  charge.  The 
index  //  (=  1,  2)  labels  the  Maxwellian.  The  extension  to 
higher-component  plasmas  would  be  done  by  simply 
including  more  terms.  The  exponentials  are  due  to  repul¬ 
sion  of  electrons  by  the  negative  potential  (Boltzmann 
factors)  and  the  square  bracketed  ion  term  represents  the 
attraction  of  positive  ions  in  the  orbit-limited  regime  of 
Mott-Smith  and  Langmuir  [1926],  which  is  a  good  ap¬ 
proximation  at  geosynchronous  altitudes.  No  approxima¬ 
tion  of  ion  mass  or  temperature  relative  to  those  of 
electrons  has  been  assumed  in  these  terms.  In  these 
expressions  it  is  assumed  that  0  <  0.  In  the  regime  0  > 
0,  the  surface  would  be  electron  attracting  and  ion  repel¬ 
ling  and  the  net  current  would  rapidly  go  to  negative 
values. 

[27]  The  symbols  in  equations  (21 )  and  (22)  are  related  to 
those  in  equations  (20a)  and  (20b)  as  follows: 

•4/ >(<5)  =  — y»jt(0)[l  -  ((S  +  (/)]exp|  -Jjr)  «  =  1.2  (24) 


=An(0)  1 


</l<? 

kTj 


n=\.2. 


(25) 


We  ignore  the  ions  for  0  near  0,  because  je  „( 0)  y,>(0) 

(n  =  1,  2).  The  flux  inequality  is  usually  valid  because  of 
the  difference  in  electron  and  ion  masses.  We  further 
assume  that  yV,i(0)  >  je% z(0),  in  order  to  illustrate  the  triple 
root  case.  The  sign  convention  in  equations  (20)  (23)  is 
that  the  incoming  electron  flux  is  negative. 

[28]  The  admissible  values  for  the  spacecraft  potentials 
are  given  by  the  zeros  (roots)  of  the  total  flux  (i.e.,  given  by 
the  current  balance  equation  at  steady  state) 


JT(d>)  =  0.  (26) 

Let  us  now  examine  the  nature  of  the  function  Jr(0)  by 
varying  0.  At  0  =  0,  7^0  =  0)  is  dominated  by  the 
contribution  from  the  electron  term  of  the  first  Maxwellian 
distribution 


M<t>)  =  -4j(0)[1  -  (^  +  '/>]exp(-^-) 


+y  1.2(0) 


<7«-0 

*7/2 


M0)~-jeA0)[\  -05  +  7/)].  (27) 

Sinee  T]  <  T*  (equation  (15)),  there  are  more  outgoing 
electrons  than  incoming  ones,  thereby  rendering  {6  +  7;) 
greater  than  unity  and  therefore  Jj{ 0)  is  positive.  This 
situation  is  depicted  in  Figure  3  (top). 
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Figure  3.  A  (top)  triple-root  situation  and  (bottom)  triple- 
root  jump  in  spacecraft  potential  0  for  a  double-Maxwellian 
distribution  of  electrons  in  spaee.  The  roots  arc  at  the  zeros 
of  the  total  flux.  The  Jr(0)  curve  for  0  >  0  follows  a 
different  equation  for  electron  attraction  and  ion  repulsion. 


[2i>]  As  the  magnitude  of  0  (<0)  increases,  the  exponen¬ 
tials  in  equations  (21)  and  (22)  beeome  increasingly  influ¬ 
ential.  Since  Tx  <  T2 ,  the  exponential  in  equation  (21) 
decreases  faster  than  that  in  equation  (22)  and  <  0) 
decreases.  Eventually,  at  sufficiently  large  magnitude  of  0 
(<0),  <  0)  is  dominated  by  the  contribution  from  the 

electron  term  of  the  second  Maxwellian  distribution.  Since 
7\  >  T*,  the  {<*>  +  7/)  term  in  equation  (27)  is  less  than  unity, 
rendering  J7{0  <  0)  negative  (Figure  3  (top)). 

JT(o  <  0)  =  -y,.2(0)(l  -  {<$  +  7/)]  exp  •  (28) 

As  the  magnitude  of  o  (<0)  increases  further,  eventually 
the  ion  terms  in  equations  (21)  and  (22)  will  beeome 
increasingly  important.  This  is  because  more  and  more 
positive  ions  are  attracted  by  the  increasingly  negative 
potential  0(<O).  This  results  in  the  total  fluxes  J7(0  <  0) 
rising  to  positive  values  (Figure  3)  and  subsequently. 

[30]  As  the  environment  parameters  change  in  time,  the 
flux-voltage  curve  adjusts  accordingly.  If  the  net  flux 
decreases  or  increases  through  zero,  a  pair  of  neighboring 
roots  ean  sometimes  eoalesee  and  disappear.  If  this  happens, 
a  jump  in  spacecraft  potential  to  the  remaining  root  would 
occur  (Figure  3  (bottom)).  In  a  general  plasma,  there  may  be 
more  than  three  roots.  The  even  roots  are  stable  and  the  odd 
roots  are  unstable  [Lai,  1991a]  If  the  dynamic  plasma 
changes  such  that  two  roots  eoalesee,  the  jump  is  between 
two  remaining  stable  roots. 

[3 1  ]  The  prerequisites  for  a  triple  root  jump  are:  a 
positive  flux  near  the  origin  (0  <  but  ^0),  with  a  decrease 
to  negative  flux  at  middle  potentials  0(<O),  followed  by 
the  rise  to  positive  flux  in  the  limit  (0  <$:  0).  If  the  positive 
flux  at  low  potentials  0(<O)  goes  negative  or  the  negative 
flux  at  the  middle  potentials  goes  positive,  two  of  the  roots 


ean  disappear  and  a  jump  to  the  remaining  stable  root 
tnggered. 

6.  Double-Maxwellian  Distribution  With 
a  Low-Temperature  Component 

[32]  Now  let  us  look  at  spacecraft  charging  when  there 
is  a  double-Maxwellian  plasma  with  a  component  that 
has  temperature  less  than  the  antieritieal  level.  For 
brevity,  we  label  the  electron  terms  by  their  uncharged 
strength  jn  =  |JC,„(0)|  where  n  =  1,  2  corresponding  to 
the  wth  temperature. 

6.1.  Case  1 

[33]  Consider  a  double-Maxwellian  distribution  in  which 
the  two  temperatures  To  and  T,  satisfy  the  following 
inequality: 

Tq  <  T4  <  T\  <  T*  (29) 

We  initially  assume  that  at  0  =  0,  the  dominant  contribution 
is  from  /o,  whose  temperature  T()  satisfies  equation  (29),  and 
therefore  Jj{ 0)  is  negative.  As  the  magnitude  of  0(<O) 
increases,  the  exponential  term  multiplying  y„  decreases 
faster  than  that  multiplying  jx.  Eventually,  the  seeond 
distribution  (T  =  T0  gains  the  upper  hand.  The  seeond 
distribution  has  TA  <  T\  <  T*  (equation  (29))  and  therefore 
J t^0  <  0)  ean  be  pushed  from  negative  to  positive.  As  the 
magnitude  of  0(<O)  increases  further,  eventually  the 
incoming  ion  flux  will  win  beeause  a  high  negative  voltage 
attracts  abundant  positive  10ns.  Therefore,  J7<0  <C  0)  goes 
more  positive.  We  eonelude  that  the  situation  of  equation 
(29)  does  not  feature  a  triple  root.  There  is  only  a  single 
negative  root,  which  arises  when  the y0  term  pushes  the  net 
flux  down  to  negative  values.  This  situation  is  illustrated  in 
Figure  4f.  If  the  jo  term  is  not  strong  enough  to  drive  the  net 
flux  to  negative  values,  the  root  would  form  at  positive 
potential.  If  the  j\  term  is  negligible,  the  flux  curve  would 
match  the  ion  line  until  the  j0  term  eomes  in.  On  the  other 
hand,  if  at  0  =  0,  the  dominant  contribution  is  from  j\%  then 
the  current  is  initially  positive  and  stays  positive.  There  is 
again  no  triple  root  situation  and  a  single  root  would  form  at 
positive  potentials. 

[34]  Consider  now  the  ease  where  the  ambient  plasma 
changes  such  that  the  j\  flux  increases.  The  net  flux  is  lifted 
higher  and  any  existing  negative  root  would  move  toward 
the  origin  (less  negative). 

6.2.  Case  2 

[35]  Consider  a  double-Maxwellian  distribution  in  whieh 
the  two  temperatures  T0  and  T2  satisfy  the  following 
inequality: 

T()  <  Ta  <  T*  <  T2.  (30) 

We  initially  assume  the  first  Maxwellian  distribution  be 
more  abundant.  At  0  =  0,  the  dominant  contribution  is  from  0 
jo,  whose  temperature  T0  satisfies  equation  (30)  and 
therefore  Jj{ 0)  is  negative.  As  the  magnitude  of  0(<O) 
increases,  the  first  exponential  (equation  (21))  loses  its 
influence  compared  with  that  of  the  seeond  (equation  (22)). 
Eventually,  the  seeond  Maxwellian  gains  the  upper  hand. 
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Figure  4.  (a)  Current-voltage  behavior  in  a  single-Maxwellian  plasma.  The  inputs  are  listed  in  Table  4 
(ease  a),  (b)  Current- voltage  behavior  in  a  single-Maxwellian  plasma.  The  inputs  are  listed  in 
Table  4  (case  b).  (c)  Current-voltage  behavior  in  a  single-Maxwellian  plasma.  The  inputs  are  listed 
in  Table  4  (ease  e).  (d  g)  Current- voltage  behavior  in  double-  or  tnple-Maxwellian  plasma.  The 
inputs  are  listed  in  Table  4  (eases  d  g). 
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However,  the  second  Maxwellian  distribution  (T2  >  T*)  also 
favors  negative  flux.  Therefore,  both  Maxwellian  distribu¬ 
tions  behave  in  the  same  way:  both  jo  and  j2  give  negative 
flux  and  there  is  no  competition  between  the  distributions. 
As  d>  (<0)  increases  further,  the  incoming  ion  flux 
dominates  because  high  negative  voltages  attract  positive 
ions.  We  conclude  that  in  this  ease  there  is  no  triple-root 
situation.  There  is  only  a  single  negative  root,  arising  when 
the  term  jo  plus  the  term  j2  push  the  net  flux  to  negative 
values.  The  characteristic  curve  is  shown  in  Figure  4e.  If 
these  combined  terms  are  not  strong  enough  to  drive  the  net 
flux  to  negative  values,  the  root  would  form  at  positive 
potential. 

[36]  If  j2  increases,  as  the  plasma  evolves,  the  net  flux 
decreases.  A  root  at  negative  potential  would  move  toward 
more  negative  potentials.  If  j2  is  strong  enough  to  produce  a 
negative  root  without  a y0  contribution,  the  root  can  move  to 
high  negative  potentials.  This  is  the  normal  ease  for  charg¬ 
ing  above  the  critical  temperature  (Figure  4a). 

7.  Triple-Maxwellian  Distribution  With  a 
Low-Temperature  Component 

[37]  Consider  a  triple-Maxwellian  distribution  in  whieh 
the  three  temperatures  satisfy  the  following  inequality: 

To  <  Ta  <  h  <  T*  <  7V  (31) 


The  equation  (19)  is  now  generalized  to 

(32) 

where  the  terms  have  the  forms  given  in  equations  (21) 
(23).  As  in  section  6,  we  assume  that  each  of  the  electron 
fluxes  is  about  2  orders  of  magnitude  higher  than  the 
corresponding  ion  fluxes.  Let  us  assume  further  that  the  first 
Maxwellian  electron  flux yeO(0)  is  stronger  than  the  seeond 
and  third  electron  fluxes.  Accordingly,  at  0  =  0,  the  total 
flux  Jr(0)  is  negative,  and  a  representative  point.  A,  is 
shown  in  Figure  5. 

[38]  As  <£(<0)  increases  in  magnitude,  the  low- 
temperatures  Tel}  in  the  denominator  of  the  exponential 
reduces  the  low-energy  flux  jeja(<t>)  faster  than  the  middle 
energy  flux  je%\(<p).  Eventually,  jet\(<j>)  exceeds  je.M)  and 
can  dominate.  Since  the  temperature  Te  \  of  je  i(0)  is  below 
the  critical  temperature  F*,  the  contribution  is  positive.  In 
other  words,  the  total  flux  Jj{(p)  can  be  driven  positive.  A 
representative  point  B  is  shown  in  Figure  5. 

[39]  As  <p(<0)  increases  further  in  magnitude,  the  expo¬ 

nential  of  in  equation  (32)  will  diminish  faster  than 
that  of  J2{4>)  in  equation  (32)  because  Te  X  <  Tc  2.  When 
J2(<P)  dominates,  charging  to  negative  voltages  will  oeeur 
because  T2  exceeds  the  critical  temperature  T*.  As  a  result, 
the  total  flux,  can  go  negative.  A  representative  point 

C  is  shown  in  Figure  5. 

[40]  As  d>(<0)  increases  even  further  in  magnitude,  even¬ 
tually  the  increasing  highly  negative  potential  0  will  attract 
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Figure  5.  A  triple-root  situation  for  a  tnple-Maxwellian 
distribution  of  electrons  in  space.  As  the  space  environment 
changes,  the  distributions  change.  Suppose  the  spacecraft 
surface  potential  is  at  the  first  root  0O  on  the  right  (near  0) 
initially.  Let  the  point  near  B  change  in  such  a  way  that  it 
sinks  as  time  evolves.  When  it  sinks  below  =  0,  the 
first  two  roots  0O  and  coalesce  and  disappear  together, 
leaving  only  the  third  root  02  at  the  far  left  (high  0  <  0).  The 
spacecraft  potential  would  jump  to  the  third  root  02. 

so  many  incoming  ions  that  the  ion  flux  will  dominate  over 
the  fluxes  of  electrons.  When  the  incoming  ion  flux  dom¬ 
inates,  the  total  flux  Jjicfi)  will  be  positive.  Further  increases 
in  the  magnitude  of  0  (<0)  will  not  alter  the  dominance  of 
the  ions.  A  representative  point  D  is  shown  in  Figure  5. 

[41]  The  curve  in  Figure  5  illustrates  a  possible  initial 
flux-voltage  characteristic  compatible  with  the  plasma  tem¬ 
perature  condition  (equation  (3 1 )).  As  the  plasma  evolves  in 
time,  the  flux  jx  could  decrease  ory2  could  increase,  causing 
the  positive  bump  in  the  flux-voltage  curve  (at  point  B)  to 
drop  down  to  negative  current.  If  the  surface  was  originally 
lying  at  the  root  between  B  and  A,  it  would  jump  to  the  high 
negative  root  between  points  C  and  D, 

[42]  In  Figure  6  we  show  a  ease  where  the yj//2  ratio  is  low 
enough  to  lower  the  net  flux  near  the  origin  to  slightly 
negative  values.  It  is  evident  that  no  triple  root  is  possible, 
but  that  one  could  occur  if  there  was  no  /0  component  pulling 
down  the  curve  at  the  origin  (so  the  plot  looks  like  Figure  3 
(top)).  It  is  in  special  cases  like  this  one  that  the  T0  term  makes 
a  difference  in  the  possible  appearance  of  triple  roots. 

[43]  This  discussion  illustrates  that,  depending  on  the 
relative  magnitudes  of  the  fluxes,  a  triple  root  situation 
can  develop  when  the  plasma  temperatures  span  the  critical 
and  antieritieal  temperatures.  The  antieritieal  component 
produces  a  negative  down  turn  in  the  flux-voltage  charac¬ 
teristic  curve,  near  the  origin.  It  only  slightly  modifies  the 
triple  root  situation  relative  to  the  double  Maxwellian  ease 
and  the  formation  of  triple  roots  still  depends  primarily  on 
the  relative  strengths  of  the  jx  and  j2  fluxes.  A  more 
quantitative  numerical  treatment  of  the  threshold  for  the 
triple  root  situations  is  described  in  Appendix  A. 

8.  Antieritieal  Temperature  and  Low-Voltage 
Charging 

[44]  In  this  section  we  describe  the  possible  shape  of  a 
low-voltage  surface  charging  curve  plotted  against  the 


average  electron  temperature.  We  will  assume  that  the 
spacecraft  is  differentially  charged  and  that  there  is  a  cloud 
of  low-energy  electrons  present.  The  electron  cloud  could 
form  because  of  secondary  electrons  trapped  locally  or 
because  of  returning  photoelectrons  in  sunlight.  The  mech¬ 
anisms  of  low-energy  electron  cloud  formation  are  outside 
the  scope  of  this  paper. 

[45]  We  consider  a  plasma  distribution  with  three  possible 
Maxwellian  components,  with  relative  flux  strengths  J0,  Ju 
J2.  The  charging  is  discussed  with  respect  to  the  average 
electron  temperature.  Let  the  overall  plasma  temperature  T 
be  given  by  the  flux-weighted  expression 

j  _  jo Fq  +ji  7j  +72  72  ^3) 

jo  +/ 1  +72 

Let  us  look  first,  at  the  lowest  average  temperature  region. 
The  flux  is  dominated  by  electrons  with  temperature  in  the 
antieritieal  region.  The  potential  is  set  mainly  by  a  balance 
between  the  low-temperature  electrons  with  strength  J()  and 
the  ions.  This  situation  is  described  in  Figure  4c,  with 
negligible  Jx  contribution.  The  electrons  add  a  negative 
contribution  near  the  origin  and  if  the  strength  J0  is 
sufficient  to  turn  the  net  flux  into  negative  values,  there  will 
be  a  root  formed  If  the  electron  flux  is  not  strong  enough, 
the  flux  will  be  positive  at  the  origin,  and  the  root  will  form 
at  positive  voltage.  Let  us  now  consider  the  ease  where  the 
average  electron  temperature  is  above  the  antieritieal 
temperature,  but  below  the  critical  temperature.  We  can 
represent  this  situation  by  adding  a  substantial  positive 
contribution  of  strength  Jx  from  a  Maxwellian  term  with 
temperature  above  antieritieal.  The  effect  of  this  term  is  to 
bump  the  JNET  curve  up  until  it  reaches  the  point  at  which  the 
low-energy  electrons  come  in  and  again  produce  a  drop  to 
negative  flux  near  the  origin  (sec  Figure  4f).  We  note  that  a 
rising  bump  due  to  Jx  would  cause  the  root  to  move  slightly 
toward  the  origin.  Consider  now  the  case  with  the  average 
electron  temperature  well  above  the  critical  temperature.  A 
dominant  electron  flux  of  strength  J2  would  pull  the  net  flux 
curve  to  negative  values.  If  there  is  any  small  low-energy  J0 
contribution  still  present,  it  would  have  a  little  effect  on  the 
root.  The  transition  from  the  situation  of  Figure  4f  to  the 
present  case  (of  rising  electron  temperature  resulting  in  J2 
dominating),  is  shown  in  Figure  6.  For  the  purpose  of 
illustrating  the  transition  from  Figure  6  (top)  to  Figure  6 
(bottom),  we  have  highlighted  the  JNET  curve.  There  would 
be  an  approximately  linear  rise  in  the  charging  potential 
above  the  critical  temperature  as  can  be  predicted  from  a 
simple  single  Maxwellian  model  [ Lai  and  Tautz ,  2006]. 

[46]  The  ease  where  the  average  electron  temperature  is 
approximately  equal  to  the  critical  temperature  can  be  quite 
complex.  In  this  region  one  expects  to  get  both  J\  and  J2 
contributions  and  the  result  is  sensitive  to  the  ratio  J\IJ2.  If 
J\  is  stronger  than  J2,  a  triple  root  situation  can  occur,  as 
depicted  in  Figure  5.  If  Jx  is  weaker  than  J2  one  can  obtain 
the  curve  shown  in  Figure  6  (bottom),  where  the  Jx  term  is 
not  strong  enough  to  push  the  net  flux  to  positive  values  and 
a  single  root  forms.  In  this  complex  ease,  the  threshold  for 
the  triple  root  formation  depends  on  the  temperatures  as 
well  as  the  ratio  J\(J2  and  is  not  very  sensitive  to y0*  1°  the 
limit  T\  T2  =  T*,  the  high-temperature  electron  fluxes 
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Figure  6.  Transition  from  (top)  a  triple-root  situation  to  (bottom)  a  single-root  situation.  The  inputs  arc 
AT0  =  0.1 1,  kT2  =  0.3,  and  AF3  =  0. 1  keV.  The  density  inputs  are  //<>  =  4,  nx  =  0.4  (in  Figure  6  (top)),  ny  = 
0.2  (in  Figure  6  (bottom)),  and  n2  -  0.1  em  \  The  ion  temperature  kTt  =  1  keVand  the  ratio  of  ion  eurrent 
to  electron  eurrent  is  taken  as  0.01. 


tend  to  eaneel  and  one  reverts  to  the  previous  low-energy 
ease  (Figure  4a). 

[47]  For  comparison,  we  have  grouped  together  the  com¬ 
plete  set  of  seven  different  eases  of  eurrent- voltage  behav¬ 
iors  in  Figures  4a  4g.  The  equations  used  are  equations 
(28)  and  (29)  or  their  respective  simplifications  for  double 
and  single  Maxwellian  eases.  The  density  and  temperature 
inputs  are  listed  in  Table  4. 

[4x]  We  remark  that  near  F*,  the  eurrent- voltage  curve 
tends  to  be  nearly  flat.  A  slight  “pulling  up”  and  a  “push 
down”  near  the  zero  crossing  may  make  a  triple  root 
situation  to  appear.  For  a  triple  root,  we  need  a  positive 
flux  J\  at  T]  <  T*  in  comparison  with  a  negative  flux  J2  at 
T2  >  7*| .  One  can  write  7*1  =  7*  Ai  and  T2  -  7*  +  A2. 
Sinee  A|  and  A2  do  not  have  to  be  equal  or  small,  it  is 
difficult  to  make  a  general  statement  on  the  triple  root 
location.  Suppose  a  triple  root  situation  forms  as  a  result. 


The  corresponding  signed  fluxes  J{)  and  J2  are  negative 
beeause  of  our  temperature  assumptions.  The  flux  j0  can  be 
smaller  than  j\  and  j2,  beeause  its  velocity  is  much  slower 
than  the  others.  In  the  average  temperature  F,  the  flux  j{)  ~  0 
in  equation  (33),  so  that 

rJiUlMl.  (34) 

j  I  +72 

If  j\  &  72,  equation  (34)  becomes 

T=i-(rt  +  r2),  (35) 

and  assuming  A|  ~  A  2  or  small  with  respeet  to  F*,  then 
equation  (35)  becomes  F  %  7*. 

[49]  In  Figure  7  we  show  empirical  data  for  the  negative 
of  the  measured  spacecraft  potential  versus  average  electron 
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Table  4. 

Electron  Density  Input  in 

Case  Studies3 

Case 

«o 

"1 

«2 

a 

0.0 

0.0 

0.1 

b 

0.0 

0.1 

0.0 

c 

1.0 

0.0 

0.0 

d 

0.0 

0.3 

0.1 

c 

4.0 

0.0 

0.1 

f 

4.0 

0.2 

0.0 

£ 

4.0 

0.275 

0.1 

“Input  densities  «o,  nj  and  n2  (cm'  3)  for  three  Maxwellian  distributions  with 
temperatures  T0  <  T4  <  T\  <  7*  <  72,  respectively.  For  the  calculations  in 
Figures  4a  4g,  we  have  used/: 70  0.011  keV,  kTA  =  0.014  kcV,A7i  =0.3  keV, 
kT*  0.5  keV,  and  kT2  =  0.9  keVs.  The  kapton  values  of  TA ,  7*,  secondary 
electron  coefficient,  and  backscattcrcd  electron  coefficient  arc  used. 

temperature  of  the  ambient  plasma,  obtained  on  a  Los 
Alamos  National  Laboratory  (LANL)  satellite  under  eclipse 
conditions  [Thomsen  et  ai ,  1999].  We  have  provided  above 
a  qualitative  physical  interpretation  to  the  data.  Only  two 
electron  distributions  plus  ions  are  needed  to  describe  the 
main  features  of  Figure  7.  At  low  average  temperatures,  low 
negative  charging  ean  oeeur.  As  the  average  temperature 
rises  above  the  anticritical  level  (here  estimated  at  —  10  eV), 
positive  flux  comes  in  and  lifts  the  net  flux  curve,  shifting 
negative  roots  to  smaller  values.  This  shift  to  smaller 
negative  charging  voltages  is  seen  statistically  in  Figure  7. 
As  the  average  temperature  goes  above  the  entical  temper¬ 
ature  (~300  eV),  the  negative  flux  begins  to  dominate  and  a 
linear  response  to  high  (negative)  potentials  is  seen.  Triple 
root  configurations  can  occur  when  the  average  temperature 
is  near  the  critical  temperature  and  the  high  (negative) 
potentials  seen  in  Figure  7  at  about  300  eV  have  the 
signature  of  triple  root  jumps. 

9.  Conclusion 

[50]  For  many  surface  materials  an  anticritical  tempera¬ 
ture  exists  along  with  the  well-known  entieal  temperature 


for  onset  of  negative  charging.  The  anti  entical  temperature 
allows  low-voltage  negative  charging  to  oeeur  below  it  and 
suppresses  the  charging  above  it.  These  two  temperatures 
oeeur  because  of  the  characteristic  shapes  of  the  secondary 
and  baekseatter  yield  eurves  and  the  theory  for  eaeh  is 
similar  (eurrent  balanee).  While  the  entieal  temperatures  for 
materials  charging  are  of  the  order  of  a  kilovolt,  the 
antieritieal  temperatures,  if  they  exist,  are  typieally  less  than 
a  few  hundred  volts.  There  is  a  large  uncertainty  in  the 
antieritieal  temperature  due  to  the  unknown  behavior  of  the 
electron  baekseatter  as  the  incident  energy  goes  to  zero. 
Reliable  values  of  antieritieal  temperature  await  the  devel¬ 
opment  of  more  accurate  electron  baekseatter  measure¬ 
ments.  In  theory,  a  low-energy  enhancement  in  the 
baekseatter  yield  moves  the  anticritical  temperature  to 
smaller  values  and  it  vanishes  if  the  baekseatter  is  too 
strong.  Due  to  the  low  energies  involved,  approximate 
values  of  the  antieritieal  temperature  ean  be  obtained  by 
Taylor  expansions,  but  these  are  also  subject  to  the  same 
baekseatter  uncertainties.  We  have  looked  systematically  at 
low-energy  charging  with  a  double  and  triple  Maxwellian 
plasma.  The  low-energy  charging  versus  cleetron  tempera¬ 
ture  is  desenbed  qualitatively  and  is  in  approximate  agree¬ 
ment  with  a  sample  of  LANL  measured  data.  Triple  root 
situations  are  also  diseussed  and  a  rough  signature  is  given; 
anomalous  charging  oeeurs  just  below  the  entieal  temper¬ 
ature.  Charging  points  having  this  approximate  tnple  root 
signature  are  also  seen  in  the  sample  data. 

Appendix  A:  Thresholds  for  Triple-Root 
Formation 

[5 1 ]  We  consider  a  triple  Maxwellian  plasma  with  tem¬ 
peratures  T0,  7j,  and  T2  which  satisfy 

Tq  <  rA  <  T\  <  r  <  r2l 


1994-084  January  1997 


Average  Te 

Figure  7.  LANL  spacecraft  charging  versus  average  electron  temperature.  The  data  is  for  eelipse.  A 
rough  estimate  of  the  anticritical  temperature  k  TA  is  10  eV,  and  a  rough  estimate  of  the  critical  temperature 
kT*  is  300  eV.  (Courtesy  of  M.  F.  Thomsen.) 
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where  TA  and  T*  are  the  anticntieal  and  critical  temperatures 
for  charging.  We  can  write  the  total  incoming  flux,  in  the 
orbit-limited  regime,  as 


yr(c»=yoexp(-g)+y,cxp( 

+  V2exp(-^)+7,('-^),  (Al) 


For  the  purposes  of  numerical  analysis,  wc  can  regard  these 
conditions  as  two  equations  in  the  unknowns  f0  =  j\y/j\  and 
fi  =  jitfu  with  To*  7i,  T2 ,  qe0 1,  qe<j>2  as  fixed  parameters.  A 
simple  matrix  inversion  then  yields  the  solutions 


JL(u\  i  (  cxp(— <?i*c!>|/AT, )  \ 
kT\  \  c\p(—qe02/kT\ ) ) 


(A5) 


where  (A/)  is  the  2  x  2  matrix 


where  qe  is  the  electron  charge,  q,  the  ion  charge,  and  0  <  0 
the  surface  potential.  The  positive  ion  flux  Jt  given  in  the 
last  term  of  equation  (Al )  is  typically  much  smaller  than  the 
electron  fluxes,  at  low  charging  potentials,  and  we  shall 
neglect  it  in  the  (necessary)  threshold  condition.  The  JN 
(N  -  0,  1,  2)  are  the  signed  electron  fluxes  at  zero 
potential.  In  terms  of  the  variables  in  the  main  text,  we 
have 


A  —  7fjv(0)(I  (6  +  ?/)), 

where  6  and  q  arc  the  corresponding  secondary  and 
backscatter  coefficients.  According  to  our  temperature 
assumption,  the  signs  of  J()  and  J2  are  negative  and  that 
of  J\  is  positive.  As  in  the  text,  wc  ean  introduce  the 
absolute  flux  strengths  yN  -  |/v|.  Wc  can  then  write 
equation  (Al)  as 

;(*)  = -yocxp(-g)  +y,exp(-^)  ,'2ex  p(-g£). 

(A2) 

In  equation  (A2),  the  electron  potential  energy  qc0  is 
positive  for  negative  charging  voltages.  The  signs  of  the 
terms  are  given  explicitly.  We  could  normalize  to  the 
central  region  flux,  making  j\  =  1,  since  this  term  must 
always  be  present  for  triple  root  formation.  If  there  is  no 
jo  flux  present,  a  necessary  condition  for  a  triple  root 
situation  is  simply  that  the  flux  be  positive  at  the  origin 
(0  -  0)  so  that  the  strengths  satisfy  j\  >  j2  [sec  Lai , 
1991a].  If  there  is  a  j0  component,  this  threshold 
condition  will  be  modified,  as  is  discussed  below. 

[52]  A  necessary  condition  for  a  triple  root  situation  to 
occur  is  that  the  net  flux-voltage  characteristic  goes  positive 
near  the  origin  and  then  sink  to  negative  values  at  larger 
qe0.  If  the  ions  are  neglected,  this  condition  ean  be  satisfied 
if  the  electron  curve  has  a  maximum  at  some  point  qe0\  and 
a  minimum  at  qe&2  >  qe<j>\.  Wc  can  write  these  extrema 
conditions  as 


,A4) 


This  procedure  yields  the  variables/)  and  f2  as  functions  of 
the  temperatures  and  potential  energies  Ify’j  is  known,  then 
by  substitution  in  equation  (A2)  wc  further  obtain  j(0\)  and 
j{0 2).  To  illustrate  the  numerics,  we  assume  nominal 
temperatures  (kcV)  for  the  regions:  AT0  =  0.01 1,  kT\  -  0.3, 
and  kT2  =  1 .0.  Wc  have  calculated  the  difference  qtxl0  = 
qe0 2  <fr0i  with  j(0\)  held  at  zero.  This  gives  the  d0 
dependence  at  the  (necessary)  threshold.  The  calculation 
indicates  that  d0  is  approximately  a  constant  at  the 
threshold.  Its  value  depends  on  the  chosen  temperatures. 

[53]  In  order  to  better  understand  the  onset  of  triple  root 
formation,  wc  examine  the  approximate  behavior  at  the 
threshold  condition. 


I  fqe0\  AT,,  we  can  neglect  the  first  term  in  (A. 7)  and  find 

•/2  =  exp(“^)'  (A8) 

where  we  have  defined  the  temperature  variable  T2\  >  0  as 


7%r, 


(7*2  T\) 


(A9) 


Now  let  us  look  for  another  approximate  expression  for  f2. 
From  the  matrix  solution,  wc  have 


(A10) 


If  qt.0\  kT„,  wc  can  again  neglect  the  first  term.  We  are 
left  with 


h 


(All) 


Combining  equation  (All)  with  equation  (A8)  wc  get 


qedO  ~  q(,<t>2  ~  qe0\ 


-  kT2\  log 


(A12) 
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Figure  Al.  Estimated  A 0  =  02  —  <t>\  versus  j2 ,  with y0  as 
fixed  parameter,  for  nominal  temperatures. 

This  expression  corresponds  to  the  numerical  result  of 
constant  q^lQ. 

[54]  In  a  physical  ease,  we  are  more  likely  to  know  the 
fluxes  jo,  j\,j  2  than  the  positions  of  the  extrema  <t>\  and  02. 
Using  equations  (A3)  or  (A4),  approximate  expressions  for 
0i  and  02  can  be  obtained  in  terms  of  temperatures  and 
fluxes.  An  approximation  for  0 j,  is  found  by  putting  the 
exponential  in  the  last  term  of  equation  (A3)  equal  to  one, 
whieh  yields 

<7,0,=  *r01  logg(^-£)),  (A13, 

where  To\  <  0  is  defined  in  the  same  way  as  in  equation 
(A9).  We  can  get  an  approximate  expression  for  02,  by 
neglecting  the  first  term  in  equation  (A4),  which  gives 

q,<h  =  *7-2i  log  (A  1 4) 

In  these  approximations,  the  log  argument  and  qe(f>\,  qe02 
are  positive,  leading  to  restrictions  on  the  flux  strengths  (jo/ 
To  +  jiiT2  < j\IT\  j2lj\  <  7%/7"i ).  When  a  representative  set  of 
parameters  To,  T\,  T2  and  jo,  j\,  j 2  are  known,  equations 
(A  13)  and  (A  14)  can  be  used  to  estimate  0l9  02  and  if  the 
condition 

•p 

qe02  -  ^0]  >  kT2 1  log--  (A  15) 

l  \ 

holds,  then,  to  the  level  of  approximation  considered, 
the  system  is  above  the  threshold.  If  we  again  let  the 
temperatures  have  the  nominal  values  given  above,  the 
value  of  the  right  hand  side  of  equation  (A  15)  is  found  to 
be  0.52. 

[55]  Using  the  approximate  formulae  (equations  (A  13) 
(A  15))  we  can  estimate  the  threshold  behavior  directly  as  a 
function  of  the  flux  strengths.  In  Figure  Al,  we  show  the 


variable  d0  as  a  function  of  j2,  over  three  decades  of /0  (in 
this  plot  we  assume  j\  =  1).  There  is  very  little  j0  depen¬ 
dence.  At  the  threshold  ( r/0  =  0.52)  it  is  evident  that  j2  is 
slightly  less  than  one.  Lower  values  of  j2  are  above 
threshold  and  higher  values  are  below  it.  Figure  Al  indi¬ 
cates  a  rough  signature  for  triple  root  charging:  there  should 
be  a  positive  flux  contribution  from  below  the  critical 
temperature  and  a  weaker  negative  (here  j2<  1)  flux  from 
above  the  critical  temperature.  On  a  plot  of  spacecraft 
charging  versus  the  average  electron  temperature,  the  triple 
root  potentials  would  appear  as  uncharacteristic  high  (neg¬ 
ative)  potentials  occurring  below  the  critical  temperature. 

[56]  The  threshold  relation  investigated  in  this  appendix 
is  a  necessary  condition.  To  obtain  both  necessary  and 
sufficient  conditions  one  must  include  the  ions.  Thus,  to 
get  a  triple  root  situation,  we  would  have 


j(qe<t>\)+Mqi<t>t)  >  0 

(A16a) 

+  J,(q,0i)  <  0. 

(A  16b) 

If  J,  ~  0,  we  get  back  from  equation  (Alba)  the  condition 
discussed  above  and  equation  (A  16b)  is  satisfied  identically 
since  y(02)  <  0  when  the  j2  term  is  dominant.  At  zero 
potentials  the  ion  fluxes  are  typically  2  orders  of  magnitude 
smaller  than  the  electron  fluxes,  beeause  of  the  smaller 
electron  mass,  and  the  approximation  Jt  =  0  is  a  good  one. 
We  are  here  assuming  that  we  are  below  the  asymptotic 
limit  0  — ►  oo,  where y(0)  — >  0  and  the  flux  is  dominated  by 
the  positive  ions.  A  net  flux  voltage  curve,  including  the 
ions,  is  depleted  in  Figure  5.  In  Figure  5,  point  B  is  shown 
as  the  positive  maximum  at  0|  and  point  C  is  the  negative 
minimum  at  02. 

[57]  The  above  analysis  could  be  earned  out  with  plasma 
densities  nN  instead  of  fluxes,  using  the  connecting  relation 

Jn  =  cnNy/J\(S(Tw)  +  //( 7"v )  -  1),  (A17) 

where  c  is  a  constant.  In  order  to  compute  the  secondary  and 
baekseatter  coefficients,  the  material  properties  are  needed. 
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